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A gauge invariant, non-local observable is constructed in pure gauge theory, which is identical to the gluon 
propagator in a particular gauge, permitting to define a non-perturbative parton mass for the gluon. This mass 
can be shown to be related to the IP — IS mass splitting of heavy quarkonia. Preliminary numerical results for 
3d SU(2) yield m A = 0.37(6)cf 2 , while from the lb spectrum one infers m A ~ 420 MeV for QCD. 



1. INTRODUCTION 

Knowledge of non-perturbative parton proper- 
ties is essential to understand the interplay be- 
tween high energy parton physics and low energy 
hadron physics. At finite temperature the decon- 
finement transition interpolates between the two 
regimes. The high temperature phase is expected 
to be governed by parton dynamics, which is en- 
coded in the Green functions of quarks and glu- 
ons. In general those are not gauge invariant. 

In perturbation theory, one fixes a gauge and 
studies partons directly. Although field propaga- 
tors are not physical observables, gauge invariant 
dynamical information is carried by their singu- 
larity structure. For example, pole masses defined 
from the gluon and quark propagators have been 
proved to be gauge invariant order by order in 
perturbation theory [Q. Resummation schemes 
have been designed to self-consistently compute 
the pole of the gluon propagator in three dimen- 
sions (2) , which is related to the "magnetic mass" 
of the thermal gauge theory. However, it has re- 
mained an open question whether a pole exists 
non-perturbatively. 

On the other hand, numerical gauge fixing on 
the lattice is problematic: It is difficult to avoid 
Gribov copies and fix a gauge uniquely. Most 
complete gauge fixings (e.g. the Landau gauge) 
violate the positivity of the transfer matrix, thus 
obstructing a quantum mechanical interpretation 
of the results. For quark mass computations, 
gauge fixing can be avoided by non-perturbative 



renormalization techniques. However, those rely 
on hadronic quantities through PCAC relations 
and are not applicable to gluons. A way around 
these problems has been suggested Q by em- 
ploying spatially non-local, composite operators. 
Without loss of generality, SU(2) pure gauge the- 
ory will be considered in the following. 

2. A NON-LOCAL GLUON OPERATOR 

A gauge invariant gluon operator can be de- 
fined when a complex N-plet transforming in the 
fundamental representation is available. One pos- 
sibility is to take the eigenfunctions of the spatial 
covariant Laplacian, which is a hermitian opera- 
tor with a positive spectrum, 

-{D 2 AU]) ap ff\x)=\ n f^\x), A">0. (1) 

They provide a unique mapping U — > f[U] ex- 
cept when eigenvalues are degenerate or |/| = 0. 
In simulations the probability of generating such 
configurations is essentially zero. These proper- 
ties have been used previously for gauge fixing 
without Gribov copies and to construct block- 
spins for the derivation of effective theories . 

The lowest eigenvectors are used to construct 
the matrix Q(x) = jjf^ / (2) 0)), which 

transforms as £l 9 (x) = g(x)&l(x)h* (t). We can 
now define composite link and gluon fields 

V„(x) = ^(x)U^x)fl(x + A), (2) 

AM = ^-(V,(x)-V}(x)), (3) 
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Figure 1. The W-boson mass in a 3d Higgs 
phase, computed from the standard operator 
V[(j>, U] and the non-local V[U], Eqs. @. 
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Figure 2. Lowest states using the composite link 
V[C] and Polyakov loops in YM-theory at (3 = 9 
with spatial lattice size L. 



both transforming as Of(x) = h(t)Oi(x)h' (t), 
whereas Vg(x) = h(t)V (x)rf(t + 1)._ The 
zero momentum projected time links Vb(i) = 
X) x Vb(x, t)/\ J2 X Vb(x, i) are multiplied to 
"strings" Vb(ti,i2) connecting two timeslices. 
These ingredients can be combined to the gauge 
invariant operator 

0[U] = Tr [^(x,0)Vb(0,i)A l (x,t)V/ o t (0,i)] , (4) 

which in the particular gauge Vo(t) = 1 reduces to 
the gluon propagator. In 0] the transfer matrix 
formalism was used to show that Eq. (Q) falls off 
exponentially as 
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The eigenvalues E n and eigenvectors |n L0 ) are 
those of the Hamiltonian in Laplacian temporal 
gauge, Vq(x) = 1. It has been proved ||] that 
this Hamiltonian has the same spectrum as the 
Kogut-Susskind Hamiltonian, which is obtained 
by quantizing the theory in temporal gauge 

Note that this does not imply an asymptotically 
free colored state. The temporal links in the op- 
erator indicate the presence of static charges, in 
analogy to the Wilson loop. Here the zero mo- 
mentum projection switches off the self-energy of 
the sources, which thus remain classical fields and 
do not affect the energy eigenvalues. 



First numerical tests are performed in the 
3d theory, which displays the same confinement 
properties but is much more easily accessible nu- 
merically. It is expedient to first test this new 
operator in a Higgs model, which in its broken 
phase has physical states with the quantum num- 
bers of the gluon, the W-bosons, with detailed 
results available Q. Fig. [j] compares the W- 
mass measured in these works by the standard 
operator V[U,</>] = Tr (T a ^ {x)U(x)(/)(x + A)), 
with results obtained from the composite links 
V[E/]. Full agreement is observed for different lat- 
tice spacings, which also implies that the energies 
Eq. (Q) do have a continuum limit. On the other 
hand, when simulated in the pure gauge theory, 
the corresponding masses grow linearly with the 
box size, as shown in Fig. |2|. This effect stems 
from the non-locality of the gluon operator, which 
depends on all link variables in a t-slice. On a pe- 
riodic torus and in a confining regime, it will thus 
project predominantly on torelonic states. This 
is evinced by correlations of the tracclcss part of 
spatial Polyakov loops L[V] constructed from the 
composite links, which produce identical results. 
The same observation is made for the ++ opera- 
tor Tr Vi (x) , which almost exclusively couples to 
the torelon, even when the latter becomes heav- 
ier than the lightest scalar glueball. In a confining 
dynamics, Eq. (Q) thus appears to only pick up 
torelonic states and is inconclusive regarding the 
mass scale associated with the gluon. 
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3. THE GLUON PROPAGATOR AND 
STATIC MESONS 

It is then necessary to construct an alternative 
operator coupling to the states we are interested 
in. This can be achieved by a limit procedure. 
Adding the scalar action 

S 4> [U 1 cj ) ]=Y J {-\D^(x)\ 2 +<\^)\ 2 } (6) 

X 

to the pure gauge theory results in QCD with 
scalar quarks. In the limit tuq — > oo the scalars 
become static sources propagating in time only, 
their propagator being known exactly to consist 
of temporal Wilson lines. 

Scalar and vector mesons are described by 
S(x) = ^(x)<f>(x),V(x) = lm(4,t(x)Di(x)<j>(x)), 
respectively. In the static limit the scalar fields 
do not contribute to angular momentum, nor 
can they be excited into higher quantum states 
since they are quenched. Consequently the mass 
difference rriA = My — Ms is a pure gauge 
quantity in that limit, characterizing an excita- 
tion with the quantum numbers of the gluon. 
Moreover, integrating the scalars out analyti- 
cally, one obtains for the ratio of correlators 
(V(x)V(y)) c /(S(x)S(y)) c ~ 
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In temporal gauge this reduces to the gluon prop- 
agator again, which then relates the mass differ- 
ence My — Ms to the pole mass of the gluon. 

Numerical results, again for 2+1 dimensional 
SU(2), are shown in Fig. ||. With increasing 
scalar mass, the measured glueball states Mq 
attain their pure gauge values indicated by the 
dashed lines. The scalar bound states move out 
of the spectrum, with My — Ms approximately 
constant. At the largest scalar mass one finds 
m A = 0.37(6)g 2 , or M G /m A w 4.2. 

These considerations are readily extended 
to QCD and heavy quark physics. Averag- 
ing over quark spins, one finds the splittings 
Aip_i S (cc) = 418.5 MeV, Ai P _i S (66) = 416 
MeV to be nearly equal, despite a factor of 
three difference in the charm and bottom quark 
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Figure 3. The lowest states in 3d scalar QCD. 
Mq denotes scalar glueballs, Ms.v scalar and vec- 
tor mesons. 



mass. According to the above, this splitting is 
determined by the QCD gluon pole mass, up to 
finite quark mass corrections. With the lightest 
pure SU(3) glueball at ~ 1700 MeV, one has in- 
deed Mc/rriA ~ 4.1, in striking analogy to the 3d 
SU(2) model. A calculation of the static limits is 
in progress. 
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